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Abstract 

The  structure  of  a  particulate  electrode  may  be  modelled  by  an  agglomerate-of-spheres,  called  the  AOS  model.  For  the  particulate 
Pb02  electrode  this  model  assumes  that  the  electronic  connections  between  the  particles  are  established  through  vase-like  necks.  The  size 
of  these  necks  determines  the  electronic  conductivity  of  the  electrode.  It  has  been  demonstrated  that  thermodynamically  there  is  a 
tendency  for  a  spontaneous  growth  of  the  neck  at  the  expense  of  the  particles.  Finally,  this  should  result  in  a  three-dimensional  network  of 
Pb02  strings.  In  the  present  paper  it  is  shown  that  the  transfer  of  material  from  sphere  to  neck  can  take  place  by  an  electrochemical  path, 
established  through  discharge  of  the  sphere  and  charge  of  the  neck.  Under  this  assumption  the  growth  rate  of  the  neck  is  calculated.  Using 
suitable  estimates  of  the  constants  involved  it  is  found  that  in  a  few  hours  the  neck  diameter  increases  to  ~  0.2  times  the  particle 
diameter,  however,  the  development  to  the  final  structure  may  take  many  days.  ©  1998  Elsevier  Science  S.A. 
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1.  Introduction 

The  agglomerate-of-spheres  (AOS)  model  was  intro¬ 
duced  by  Winsel  and  co-workers  [1-3]  as  a  three-dimen¬ 
sional  network  of  equally  sized  PbOz  spheres  connected 
by  necks.  The  present  work  concerns  primarily  the  elec¬ 
tronic  conductivity  perpendicular  to  the  surface  of  the 
electrode.  From  this  reason  only  one-dimensional  parallel 
strings  of  spheres  are  considered. 

The  structure  is  defined  by  three  interdependent  radii  of 
curvature  namely  the  sphere  radius,  rs,  the  radius  of  the 
cylinder,  which  can  be  inscribed  in  the  neck,  rn  and  the 
radius  of  the  toroidal  part  of  the  neck,  rv  see  Fig.  1. 

The  sphere  radius  can  be  estimated  from  specific  sur¬ 
face  area  data,  and  therefore  it  is  most  convenient  to  use  rs 
and  the  ratios 

rn  rt 

a  =  —  and  b  =  —  (1) 
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as  variables.  According  to  Ref.  [1]  a  is  also  the  ratio 
between  the  electronic  conductivity  of  the  particulate  mass 
and  the  bulk  conductivity. 

From  the  geometry  of  the  sphere-neck  configuration  it 
can  be  found  [1,4] 


b  = 


2(1  -  a) 


(2) 


By  thermodynamic  analysis  of  the  AOS  model  it  was 
found  [4]  that  the  transfer  of  PbOz  from  the  bulk  of  the 
sphere  to  the  neck  zone  causes  a  change  in  partial  molar 
free  energy 


1 

a 


(3) 


where  y  is  the  interfacial  tension  between  PbOz  and  the 
electrolyte  and  V0  the  molar  volume  of  PbOz. 

Within  the  limits  of  the  model  (0  <  a  <  1,  b  >  0)  it  was 
shown  that  AG  always  is  negative.  Thus,  if  a  path  for 
transfer  of  Pb02  from  the  sphere  to  the  neck  exists,  a  will 
continue  to  grow,  finally  causing  the  AOS  structure  to 
degenerate  into  a  cylindrical  string. 
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Fig.  1.  (a),  (b)  Geometry  of  the  AOS  structure. 


The  importance  of  this  development  originates  from  the 
recognition  that  the  electronic  conductivity  of  the  mass  is 
proportional  to  a  [1],  and  that  it  is  likely  that  in  some 
instances  the  conductivity  is  a  dominating  factor  for  deter¬ 
mining  the  capacity  of  the  Pb02  electrode  [5,6]. 

2.  Electrochemical  transfer 

Both  the  sphere  and  the  neck  establish  electrochemical 
exchange  with  an  electrolyte  containing  Pb2  +  ,  H+  and 
H20  combined  with  an  electron  source. 

The  process  consists  of  several  partial  reactions  [7].  At 
stationary  conditions  the  total  process  is 

Pb2+  +  2H20  =  Pb02Hx  +  (4  -  x)H+  +  (2  -  x)e~  (4) 

using  the  fact  that  Pb02  is  a  non-stoichiometric  oxide 
containing  protons  and  electrons  [8], 

The  equilibrium  potential  connected  with  Eq.  (4)  can  be 
written  as  a  standard  potential  e°  plus  terms  depending  on 
the  chemical  potentials  of  Pb2+  and  H+  in  the  solution 
and  H  in  the  oxide.  In  fact  only  two  of  these  are  indepen¬ 
dent  variables  [9]. 

Considered  separately  e°  is  different  for  the  sphere  and 
the  neck  as  the  change  in  the  G  function  by  adding 
Pb02H  r  is  different  because  of  the  difference  in  interfa¬ 
cial  energy.  As  there  is  chemical  equilibrium  between  the 
sphere  and  the  neck  [4]  it  is  clear  that  the  equilibrium 
electrode  potentials  differ  by  an  amount  given  by  AG/2F. 
(The  denominator  should  correctly  be  (2  —  x)F,  but  x  is 
negligibly  small.) 

Indicating  sphere  and  neck  by  ‘s’  and  ‘n’,  respectively, 
we  have  for  the  cell 

(Pb02  Hx)  J|H+,  Pb2  +  ,  H20||(Pb2OHx), 


e.m.f.  =  es°  -  en° 

=  -AD/2F 

(5) 

or  using  Eq.  (3) 

y  y0 

1  1  1  \ 

(6) 

e.m.f.  = - 

- 2 

2  Frs 

\  a  b  / 

If  an  electronic  connection  exists  between  sphere  and 
neck  the  non-zero  value  of  e.m.f.  will  cause  electrochemi¬ 
cal  transfer  of  Pb02  and  electrons 

[(2  -x)e”  +  Pb(02H,)]s  -  (Pb2  +  ,H+,H20)e  (7a) 

(Pb2+,H+  H20)e  — >  [Pb02Hx  +  (2  — x)e  ]n  (7b) 

(2  —  x)e~  ->  (2  —  x)es~  (7c) 

Eqs.  (7a)  and  (7b)  take  place  through  the  electrolyte  (e) 
and  Eq.  (7c)  through  neck/sphere  boundary. 

The  direction  of  Eqs.  (7a-c)  is  given  by  the  negative 
value  of  AG.  The  sphere  is  the  positive,  the  neck  the 
negative  electrode  in  the  cell.  Eqs.  (7a-c)  cause  the  sphere 
to  be  discharged  and  the  neck  to  be  charged. 

The  e.m.f.  is  the  driving  force  for  the  transfer  and  its 
magnitude  together  with  the  active  areas  determines  the 
current  associated  with  Eq.  (7b).  If  b  is  expressed  by  Eq. 
(2),  the  value  of  the  e.m.f.  can  be  calculated  as  a  function 
of  the  ratio  a,  using  V0,  y  and  rs  as  parameters. 


neck/sphers  ratio:  a 

Fig.  2.  E.m.f.  as  a  function  of  a  (Eq.  (6)),  rs  =  60  nm,  and  y  =  l  Jm 
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V0  is  well  known,  being  24.6  X  10-6  m3  mol-1.  Esti¬ 
mation  of  y  is  quite  arbitrary.  No  relevant  literature  data 
has  been  found.  In  Ref.  [1]  y  =  1  J  m  2  is  proposed.  This 
is  comparable  with  known  data  for  metals. 

Using  the  AOS  model,  rs  is  estimated  by  approximating 
the  Pb02  electrode  by  an  agglomerate  of  monodisperse 
spheres  and  equalizing  their  surface  area  to  the  known 
BET  area.  The  reason  for  using  the  BET  surface  for 
calculating  rs  is  that  the  key  point  in  the  work  is  the  use  of 
interfacial  energy,  which  is  proportional  to  the  total  sur¬ 
face.  Under  the  assumptions  of  the  model  this  surface  must 
be  converted  to  sphere  radii.  With  5  m2  g  1  as  a  typical 
value,  rs  comes  out  as  ~  60  nm. 

In  Fig.  2  e.m.f.  is  shown  as  a  function  of  a.  It  is 
noticeable  that  e.m.f.  decreases  to  rather  small  values  for 
a  >  0.5. 


3.  Transfer  current 

During  the  flow  of  the  transfer  current  the  e.m.f.  is 
balanced  by  the  transfer  overvoltages  at  the  PbO,  elec¬ 
trolyte  boundaries  in  combination  with  the  ohmic  voltage 
drop  between  the  sphere  and  the  neck  and  a  possible 
concentration  overvoltage.  As  an  approximation  (the  valid¬ 
ity  of  which  is  discussed  later)  only  the  transfer  overvolt¬ 
age  is  considered  significant.  Then 

e.m.f.  =  ??„  -  i?s  (8) 

where  17  stands  for  the  transfer  overvoltage  (pn  >  0,  ps  <  0). 
The  transfer  current  is  given  by  these  overvoltages  and  the 
electrochemical  active  areas  of  the  sphere  and  neck,  re¬ 
spectively.  Owing  to  the  magnitude  of  the  e.m.f.  the 
non-linear  Volmer-Butler  relation  between  the  current 
densities  and  the  overvoltages  must  be  used.  Indicating  the 
areas  for  neck  and  sphere,  respectively,  by  An  and  As  the 
currents  to  the  neck  and  the  sphere,  ia  and  is  are 
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where  the  exchange  current  i0  is  defined  with  respect  to 
the  electrolyte;  a  and  (3  are  the  anodic  and  cathodic 
transfer  coefficients. 

The  currents  in  the  electrolyte  at  the  neck  and  the 
sphere  are  numerically  equal,  but  of  opposite  sign 

- (10) 

By  combining  Eqs.  (8)— (10)  and  inserting  values  for 
e.m.f.,  An  and  As,  the  overvoltages  can  be  found.  Then, 
by  inserting  a  value  for  iQ  in  Eq.  (9)  the  current  transfer¬ 
ring  Pb02  from  the  sphere  to  the  neck  can  be  found  as 
function  of  a.  Finally,  from  the  relation  between  the 
volume  of  the  neck  and  a,  the  rate  of  the  increase  in  neck 
size  and  consequently  a,  can  be  found. 


4.  Experimental 

The  constants  /0,  a  and  /3  in  Eq.  (9)  were  determined 
experimentally. 

Two  commercial  positive  plates,  2  mm  thick,  were  cut 
into  5  cm  X  8  cm  surface  area.  The  active  mass  in  each 
was  29  g.  They  were  cycled  5  times  in  parallel  against  an 
oversize  negative  in  3.5  mol/dm3  H2S04.  The  cycling 
procedure  was:  discharge  with  24  A/kg  to  1.5  V  and  IUI 
charge  with  the  same  current  to  20%  overcharge.  They 
delivered  90-100  Ah/kg  Pb02.  To  ensure  the  presence  of 
PbS04  they  were  both  discharged  with  20  Ah/kg.  In  this 
way  ~  9%  of  the  total  amount  of  Pb02  was  converted  to 
PbS04. 

With  one  plate  as  reference  electrode  the  other  was 
discharged  and  charged  potentiostatically  in  the  same  acid 
against  a  negative  plate.  The  potential  range  was  +  16  mV 
and  the  geometrical  stationary  current  density  was  less 
than  1  A  m-2. 


5.  Results 

Initially  at  the  potential  step  the  current  was  larger  than 
noted  above.  This  is  due  to  the  hydrogen  insertion/extrac¬ 
tion  reaction  [7,9]  having  a  much  larger  rate  constant  than 


Fig.  3.  Anodic  and  cathodic  current  on  positive  plate  40  cm2.  X : 

experimental  values:  ( - )  Volmer-Butler  approximation  with 

(0  =  12.2X  10-3  A  and  transfer  coefficients  0.34  (anodic)  and  0.66 
(cathodic). 
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the  lead  dissolution/deposition  reaction  [7],  Within  5  to  10 
min  the  value  of  x  at  the  surface  is  adjusted  to  the  given 
overvoltage  and  the  current  acquires  a  constant  value 
representing  the  rate  of  the  stationary  Eq.  (4).  It  is  the  rate 
constant  of  this  reaction,  which  should  be  used  in  Eq.  (9). 
During  the  measurements  less  than  1%  of  the  Pb02  and 
the  acid  in  the  pores  was  used.  Due  to  the  high  electronic 
and  ionic  conductivities,  compared  with  the  low  exchange 
current  found,  it  is  estimated,  that  the  influence  of  the 
porous  structure  can  be  neglected. 

The  results  are  shown  in  Fig.  3.  A  reasonably  good  fit, 
shown  by  the  full  line,  is  obtained  with  a  =  0.34  and 
(3  =  0.66.  The  exchange  current  for  the  plate  came  out  as 
12.2  X  1 0  3  A,  corresponding  to  3.0  A  m  2  geometrical 
area. 

Using  the  BET  area  5  m2  g  1  as  an  estimate  of  the  true 
active  surface,  from  the  reasons  given  above,  renders  a 
value  of  the  exchange  current  density,  z'0,  for  Pb02  equal 
to  93.5  X  10~6  A  m-2. 

The  values,  in  particular  for  i0,  found  here  are  in 
considerable  disagreement  with  literature  data  [10-12].  i0 
is  considerably  smaller  than  the  values  cited  in  Refs. 
[10,11].  Better  agreement  is  found  with  the  value  given  in 
Ref.  [12]  for  a  comparable  plate.  This  disagreement  with 
literature  data  may  originate  from  the  large  difference  in 
the  rate  constants  between  the  hydrogen  insertion  [7]  and 
the  stationary  reaction,  which  is  not  discussed  in  the 
literature. 

One  may  question  the  validity  of  using  data  for  small 
overvoltages  for  the  e.m.f.  range,  which  is  relevant  for  the 
AOS  structure.  The  reason  for  this  was  that  in  order  to 
minimize  the  influence  of  the  porosity  of  the  electrode,  the 
current  density  should  be  small. 

As  will  appear  from  the  following,  the  time  scale  for 
the  growth  of  the  neck  is  inversely  proportional  to  i0. 
Thus,  the  choice  of  this  parameter  is  important,  but  for  the 
reasons  given  above  the  experimental  value  found  here  is 
used  in  the  following. 


Fig.  4.  Transfer  current  in  10-3  A  kg  1  Pb02  and  combined  active  area 
in  10~ 18  m2  as  a  function  of  a. 


The  solution  of  this  equation  for  r/n  for  a  given  e.m.f. 
requires  knowledge  of  the  areas  An  and  As.  These  are 
found  from  the  geometry  (Fig.  1)  by  suitable  integrations 
(see  Appendix  A) 


=  2  rrrs2 


1 

T+b 


(12) 


An  =  27rr^b\(a  +  fr)arcsin- 


— ) 
1  +  b  ) 
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From  these,  rjD  can  be  found  as  a  function  of  a  by  an 
iteration  procedure. 


6.  Calculations 


The  data  and  relations  given  so  far  allow  a  calculation 
of  how  the  ratio  a  should  increase  spontaneously  with 
time. 

The  first  step  is  calculation  of  e.m.f.  for  a  chosen  value 
of  a  according  to  Eq.  (6).  Then  combination  of  Eqs. 
(8) — (10)  give 
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Fig.  5.  a  as  a  function  of  time  0-200  min. 
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Fig.  6.  a  as  a  function  of  time  0-200  h. 


V  =  7rrs3£>j[(a  +  b )2  +  7>2]sin  00  —  b(a  +  b) 

X  (sin  60  cos  0Q  +  60)  —  b2sin  0oj  (16) 

with 


sin  60 


1 

1  +b 


;  cos  60  = 


a  +  b 
1  +b 


1 

6n  =  arcsin - 

0  1  +  b 


Integration  was  performed  by  a  stepwise  numerical 
procedure  using  a  as  independent  variable,  calculating  AV 
and  i  and  inserting  in  Eq.  (15).  The  results  are  shown  with 
a  as  a  function  of  time  in  Figs.  5-7. 


7.  Discussion 


The  transfer  current,  i  is  then  found  as  a  function  of  a 
using  Eq.  (9).  i  is  influenced  by  the  e.m.f.  (Eq.  (6))  as  well 
as  by  the  magnitude  of  the  two  areas  An  and  As.  The 
smaller  of  these  has  the  major  influence.  A  combined 
active  area  is  defined  as 

1  1  1 

+  T  (l4) 

Ar  increases  to  a  maximum  at  a  =  0.8  and  therefore 
counteracts  the  decrease  in  e.m.f. 

i  (expressed  as  mA/kg  PbO,)  and  Ar  as  functions  of  a 
are  shown  in  Fig.  4.  The  S  shape  of  i  is  due  to  the 
interaction  between  e.m.f.  and  Ar. 

The  transfer  current  indicates  the  rate  by  which  Pb02  is 
transferred  from  the  sphere  to  the  neck,  increasing  the  neck 
volume  and  diameter.  The  time  dr  to  increase  the  volume 
by  dV  is  then 

2  F 

dl=—dV  (15) 

1  *'0 

In  order  to  pass  to  a  as  variable,  the  neck  volume  is 
found  as  a  function  of  the  angle  d0  (Fig.  1(b))  and  a  and 
b_( see  Appendix  B) 


Fig.  7.  a  as  a  function  of  time  0-120  days. 


In  the  calculations  rs  has  been  assumed  constant.  How¬ 
ever,  in  the  range  of  interest  the  decrease  of  rs  is  com¬ 
pletely  negligible.  At  a  —  0.6  (after  500  h)  the  sphere 
radius  is  still  99%  of  the  initial. 

In  the  derivation  of  Eq.  (8)  the  ohmic  voltage  drop 
within  the  AOS  structure  was  neglected.  Due  to  the  low 
resistivity  of  Pb02  (1.2  X  10"6  fl  m)  and  the  short  dis¬ 
tance — taken  as  rs — the  ohmic  drop  is  orders  of  magni¬ 
tude  smaller  than  the  overvoltages. 

The  same  applies  to  a  possible  contribution  of  an  ohmic 
voltage  drop  in  the  acid  and  a  possible  overvoltage  caused 
by  a  concentration  gradient  between  sphere  and  neck. 
Using  the  solubility  of  PbS04  (2  mg  dm-3)  and  the 
diffusion  coefficient  for  Pb2+  (0.9  X  10-9  m2  s-1),  a 
diffusion  length  equal  to  rs  and  the  values  for  the 
transfer-current  density  at  the  neck,  it  can  be  calculated 
that  even  for  a  =  0. 1  the  concentration  overvoltage  is  less 
than  1 0 “ 3  times  the  transfer  overvoltage. 

The  largest  uncertainty  in  the  calculation  is  due  to  the 
values  ascribed  to  y  and  rs.  However,  this  does  hardly 
invalidate  the  conclusions  drawn  on  the  basis  of  the  AOS 
model. 

One  may  question  the  validity  of  the  model  for  the 
actual  Pb02  electrode  in  the  lead/acid  battery.  The  as¬ 
sumption  about  spheres  of  the  same  size  arranged  in  a 
one-dimensional  array  is,  of  course,  not  realistic.  Expan¬ 
sion  of  this  simplified  AOS  model  to  a  three-dimensional 
network  is  mathematically  possible  but  requires  assump¬ 
tions  about  the  coordination  number  in  the  model  [13]. 
This  expansion  would  influence  the  time  scale  of  the 
growth,  but  would  not  change  the  overall  picture  signifi¬ 
cantly.  More  questionable  would  be  the  influence  of  a 
statistical  distribution  of  the  sphere  sizes,  which  would 
complicate  the  mathematical  treatment.  This  problem  is 
worked  on  elsewhere.  1 


1  A.  Winsel,  personal  communication. 
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However,  it  is  clear  that  in  a  particulate  electrode 
consisting  of  one  component  only,  the  electrode  structure 
will  exhibit  concave  as  well  as  convex  surfaces  in  order  to 
have  electronic  conductivity  and  mechanical  coherence. 
The  AOS  model  should  be  considered  as  an  approximation 
to  this  fundamental  property. 


8.  Conclusions 

One  conclusion  of  this  work  is  that  even  if  there  is  only 
a  point  contact — of  unit  cell  dimensions — between  two 
particles  it  should  only  take  a  few  hours  for  the  neck  to 
grow  to  a  diameter  which  is  ~  20%  of  the  particle  diame¬ 
ter  (Fig.  5).  This  should  ensure  an  interparticle  conductiv¬ 
ity  sufficient  for  most  applications. 

After  this  initial  phase  the  growth  rate  is  slow.  On  a 
practical  time  scale  one  may  consider  the  range  0.3-0.5 
for  the  ratio  (a)  between  neck  and  sphere  diameter  as  a 
‘pseudostable’  state  (Fig.  6).  Afterwards  the  growth  rate 
becomes  extremely  slow  (Fig.  7). 

From  the  expression  for  e.m.f.  Eq.  (6),  it  can  be  seen 
that  particle  sizes  larger  than  used  here,  delay  the  growth. 
This  would  also  be  the  case  for  a  smaller  value  of  y. 

Tacitly  it  has  been  assumed  that  the  electrode  is  at  rest. 
However,  the  considerations  above  are  also  valid  during 
discharge  and  charge  conditions.  As  the  concentrations  at 
the  sphere  and  neck  surfaces  are  equal  (see  above)  the 
e.m.f.  is  caused  by  a  difference  in  the  standard  potentials 
(Eq.  (5)),  therefore  the  overvoltages  associated  with  the 
transfer  current  adds  to  the  overvoltage  caused  by  the 
discharge  or  charge  currents.  The  AOS  overvoltage  at  the 
neck  is  positive,  and  calculation  shows  that  in  the  range  of 
interest  it  constitutes  more  than  90%  of  the  AOS  e.m.f. 
Therefore,  in  particular  for  small  values  of  a  (Fig.  2), 
during  discharge  the  total  overvoltage  and  thus  the  current 
density  at  the  neck  is  smaller  than  the  current  density  at 
the  sphere.  At  charge  the  current  density  on  the  neck  is 
larger.  The  exponential  dependence  of  the  current  on  over¬ 
voltage  (Eq.  (9))  for  overvoltages  in  the  practical  range  of 
discharge/charge  currents  amplifies  this  effect.  From  this 
reason  the  growth  of  the  neck  should  be  favored  by 
cycling. 

The  development  of  the  electronic  resistance  of  the 
positive  mass  has  been  reported  [2,3],  During  most  part  of 
the  discharge  this  resistance  remains  approximately  con¬ 
stant,  which  agrees  with  the  reasoning  above.  However  the 
resistance  increases  abruptly  at  the  end  of  discharge.  Also 
during  repeated  cycling  the  resistance  increases  [2],  This 
cannot  be  explained  easily  by  the  AOS  theory  as  presented 
here.  Attention  is  drawn  to  the  degradation  phenomena 
(sludge  formation)  associated  with  the  expansion  of  the 
electrode  [14].  This  degradation  apparently  takes  place, 
although  the  increase  in  neck  size  should  improve  the 
coherence  of  the  electrode,  and  therefore  counteract  a 
mechanical  disruption. 


9.  List  of  symbols 

An  r  s  Active  areas  of  neck,  neck/sphere  and  sphere, 
respectively 

a  ratio  rn/ rs  (see  below) 

b  ratio  rt/rs  (see  below) 

e.m.f.  electromotive  force 

F  Faraday  constant 

AG  change  in  partial  molar  Gibbs  free  energy 

i,  zn,  zs  transfer  currents 

i0  exchange  current 

R  gas  constant 

rn  radius  of  neck  (Fig.  1) 

rs  radius  of  sphere  (Fig.  1) 

rt  radius  of  toroidal  part  of  neck 

x  stoichiometric  coefficient 

V  volume  of  neck 

V0  molar  volume  of  Pb02 

Greek  letters 


ot,  /3 
y 

V 

e,  e0 


anodic,  respectively  cathodic  transfer  coeffi¬ 
cients 

interfacial  tension 
standard  potential 
overvoltage 
angles  (Fig.  1) 


Subscripts 


n  subscript  for  neck 

s  subscript  for  sphere 


Appendix  A 


A.l.  Active  area  of  the  sphere 


With  reference  to  Fig.  1(a)  x  is  the  distance  from  the 
centerline  to  a  point  on  the  periphery  of  the  sphere,  and  ds 
is  the  angle  between  radius  to  that  point  and  the  centerline. 
Then 

dAs  =  2Trrsd0s  (Al/1) 

and 

x=rssin0s  (Al/2) 

Then 


As  =  —  2 77 rs2 jfj2 sin 0S  d0s  =  27rr2  cos6»0  (Al/3) 


where  90  is  the  angle  to  the  point  where  the  neck  touches 
the  sphere;  cos  0o  is  found  from  the  figure 


cos  0q 


G 

A  +  T 


1 

1  +b 


(Al/4) 


Then 

A*  =  2lTr'T+b  (Al/5) 
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A.2.  Active  area  of  the  neck 

With  reference  to  Fig.  1(a)  6n  is  the  angle  between  the 
basis  line  and  a  radius  in  the  toroid  to  a  point  at  the  neck. 
Then 

dAn  =  2iryrld0n  (Al/6) 

and 

y  =  rn  +  rt-rtcosda  (Al/7) 

Then 

An  =  2^>\[{rn  +  rtJo°den  ~  T/o°COS^  d0„]  (Al/8) 

where  60  is  the  complement  to  the  0O  used  above.  Intro¬ 
ducing  0Q  =  arcsin[l/(l  +  b)\  and  a  and  b  one  finally 
gets 

A„  —  2tt  r2b  (  a  +  b)arcsinl - ) - (Al/9) 

n  s  J  \l+b)  1  +b  V  '  ’ 


Appendix  B 

B.l.  Volume  of  the  neck 

With  reference  to  Fig.  1(b)  a  secant  perpendicular  to  the 
centerline  at  point  O  is  drawn  in  the  distance  y  from  the 
basis  line.  The  secant  cuts  the  sphere  at  S  and  the  neck  at 
N. 

With  OS  =  jcs  and  ON  =  xn  a  volume  element  in  the 
neck  is 

dV=  ir(x2n  ~x2)dy  (A2/1) 

Introducing  as  integration  variable  the  angle  0  between 
the  radius  in  the  toroid  to  N  and  the  basisline  we  get 

xn  =  rn  +  rt  —  rt  cos  0  =  rfa  +  b  —  b  cos  9)  (A2/2) 

and 

y  =  r,  sin0  =  rsbsin0=>  dy  —  rsb  cosd  dd  (A2/3) 


and 

x\  =  r2  —  (  rs  —  y)2  =  r2(2b  sin0  —  b 2  sin#2)  (A2/4) 
then 

—  2 b(a  +  b)  f  cos2d  dd  —  2b  f  sin  6  cos  6  dO]  (A2/5) 
Integrating  from  0  to  0o  gives 

V  =  Trr2b{  [(n  +  b)2  +  Z?2 ]  sin 0O  —  b(a  +  b) 

X(sin0o  cos0o  +  60)  —  b2smd0  (A2/6) 

with 


1  a  +  b  1 

sin  0Q  = - ;  cos  0O  = - ;  0O  =  arcsin - 
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